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Abstract 

In this paper, a new simple but yet efficient spectral expression of the frequency -limited 
'H.2-norm, denoted T^^-norm, is introduced. The proposed new formulation requires the 
computation of the system eigenvalues and eigenvectors only, and provides thus an alternative 
to the well established Gramian-based approach. The interest of this new formulation is in 
three-folds: (i) it provides a new theoretical framework for the H2, ^-norm-based optimization 
approach, such as controller synthesis, filter design and model approximation, (ii) it improves 
the "H2.^-norm computation velocity and it applicability to models of higher dimension, and 
(Hi) under some conditions, it allows to handle systems with poles on the imaginary axis. Both 
mathematical proofs and numerical illustrations are provided to assess this new "H2,u-norm 
expression. 

1 Introduction 

Norms associated to Multiple Inputs Multiple Outputs (MIMO) Linear Time Invariant (LTI) dy- 
namical models, such as rloo and - H 2 _n orms, take a predominant role in many control optimization 
problems. Indeed, these norms are often used as minimization cost for controller [51 [TTJ [3] and 
observer design, filter optimization or large-scale model approximation [51 116j . More specifically, 
due to practical considerations, frequency-limited versions of these norms - i.e. their computation 
over a limited frequency range - are very meaningfull for engineers. Indeed, (i) some frequencies 
are physically meaningless and can be regarded as uncertainties, (ii) for vibration control, some 
frequencies are more specifically of interest and (Hi) in practice, the bandwidth of actuators is 
limited making some frequencies irrelevant for control purpose. These points are specifically true 
in many applicative fields such as aerospace |13j and more particularly for gust load alleviation 
PJ 0] , very large-scale integration systems [3; , civilian engineering [7] , where attention of engineers 
and researchers focuses on a finite frequency interval only. In these restriction of the 

T^-norm over a bounded frequency range can be more appropriate. Therefore, its computation 
plays a pivotal role in many control problems, specifically when the problem dimension increases. 

This paper focusses on the T^-norm only, and more specifically, to the new definition (and 
computation) of its frequency-limited version, denoted which denotes the evaluation of the 

T^-norm from to an upper bound u> € R+. As grounded on the spectral information of an LTI 
system, the contribution of this work stands in a new spectral expression of the frequency-limited 
rl2-norm. This new surprisingly simple rt2,u formulation, which only requires the eigenvalues and 
eigenvectors computation, (i) exhibits many advantages with respect to the standard approaches 
such as applicability to large-scale models and (ii) provides as an alternative to the Gramian-based 
approach, a new framework for "H2, w optimization procedures, grounded on spectral informations 
only. 
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This paper is divided as follows : Section [2] provides preliminary results on TL2 and %2,u- 
norms definitions and computations. In Section [3j the main result of the paper, i.e. a new 
spectral expression for the frequency-limited "Hi-norm is introduced. Then, Section [4] illustrates 
the behaviour of the T-^u-norm and deals with the numerical considerations related to this new 
spectral-based formulation. Finally, Section [5] discuses the results and potential exploitation of 
this result. 

Mathematical notations are standard : the system state is denoted x £ K™ , the MIMO LTI 
dynamical system is denoted S and H(s) := C(sl n — A)^ 1 B + D is its associated transfer function. 
The T^w-norm stands for the frequency- limited %2-norm taken over the frequency interval [0,w], 
uj € R_|_. A T and A* are the transpose and conjugate transpose of the matrix A, respectively. 
The matrix logarithm and trace operators are denoted logm(.) and tr(.), respectively. A(.) holds 
for the eigenvalue operator and A^ denotes the ith eigenvalue. Let Re(z) and Im(z) be the real 
and imaginary parts of the complex number z — a + jb, where j 2 = — 1. The residual of the 
complex valued function f(z) at Aj is denoted <j>i = lim z ^\. (z — Xi)f(z). The arctangent atan(z), 
logarithm log(z) and arccotangent acot(z) of the complex variable z are complex functions defined 
in El 



2 Preliminary results on H2 and %2,u;- norms 

Let consider n u inputs, n y outputs MIMO LTI dynamical systems of order n represented as: 

x{t) = Ax(t) + Bu(t) 

y(t) = Cx\t) + Du{t) '^ 

where A £ R nxn , B € R nxn ^, C £ E"« x ", D £ E»« x »" and its associated transfer function 
H(s) = C(sl n - A)- X B + D £ C n * xn ", the "H 2 -norm of E is defined as follows. 

Definition 1 (%2-norm). Given a continuous MIMO LTI dynamical system E described as in 
([T]), the H2-norm o/E, denoted by ||E||% 3> is given as 

\\n\l 2 --=tr(±- [°° \H(ju)\ 2 dv). (2) 



2tt _ 

When E is asymptotically stable and strictly proper, then ||E||^ 2 is finite and can be computed 
as 

||E||^ 2 =tr (BQB T ) =tr(C T PC), (3) 

where V and Q are the reachability and observability Gramians solution of the following Lyapunov 
equations 

AV + VA T + BB T = , (4a) 
A T Q+QA + C T C = 0. (4b) 

Moreover, if E has simple poles A, (i = 1, . . . , n), the T^-norm can also be computed as chap. 

5] 

IISII^^trf^^C-A^, (5) 

where A,-, <j>i are the poles and residues of H(s), respectively. To the authors' knowledge, this 
residue-based formulation is mainly used in H,i optimal model reduction since it allows one to 
derive first-order optimality conditions directly in terms of interpolation conditions [5] . 

Similarly to as Definition [I] and as firstly suggested in [2] , the frequency limited T^-norm can 
be given as follows^ 

1 This norm can for instance be used for robust performance analysis [11] and in model approximation [12) . 
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Definition 2 (%2, w -norm). Given a continuous MIMO LTI dynamical system E described as in 
([I]), the frequency-limited H2 -norm o/E, denoted ||E||^ 2lJ , is £/ie restriction of its T-l2-norm over 
a bounded frequency range [0,oj], to £ 1R + , and is given as 




As rooted on the Gramian-based %2-iiorm, for a stable and strictly proper system, the resulting 
frequency-limited %2-norm can be defined as 

ll^llk. = tr (BQ^B T ) = tr (C T V^C) , (7) 
where V u and are the frequency-limited reachability and observability Gramians 7J 



1 f u 

- / T(v)BB T T*(v)dv, (8a) 



^ = 2 

Q u = ^ I T*{v)C T CT{v)dv (8b) 

Z7T 

where T(z/) = {jvI n -A)~~ . They may alternatively be obtained by solving the following two 
Lyapunov equations, 

AV U + V U A T + W c (uj) = , (9a) 
A T Q bJ + Q u A + W {u J ) = (9b) 



where 



W c (uj) = S{lo)BB t + BB t S*{lo) , (10a) 
W (u) = S*(u)C T C + C t CS(uj) (10b) 



and 



S(u) = i- / 7( 



—logm ((A + juI n )(A - jujl n ) l ) 



(11) 



Remark 1 (General frequency interval). The frequency-limited 'H.2-norm can also be defined over 
the bounded interval [wi,^], uj 1 , lu 2 € M + , < W2, as 

m\k,^ 2] = \\n\k„-\\n\k,„ 1 - (12) 



3 Main result 

The aim of this paper is to present a new expression (up to the authors knowledge) of the %2,ui- 
norm based on the spectral informations of the system, i.e. the residues and eigenvalues and of 
the transfer function H(s), instead of the Gramian-based approach. This new simple, yet efficient, 
formulation allows new perspectives from a control, observer, filter design or model reduction point 
of view and provides an alternative to the Gramian approach for the frequency limited %-norm 
computation of large-scale models. The main result of the paper is stated below. 

Theorem 1 (Spectral expression of the T^.w-norm). Given a continuous MIMO LTI dynamical 
system described as in ([!]) of degree n with simple poles \, i — l,...,n. Let <f>i £ c n " xn " 
denote the corresponding residues of H{s) at Xi, i.e. <fii — lim s ^A, (s — Xi)H(s), i = l,...,n. 
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Suppose that the purely imaginary poles A™, k = 1,.. . , rii m (0 < n{ m < n) of H(s) are such that 
to < rain (|A™ 1 |) . Then, the frequency-limited T-L2-norm can be written as 

ll^llk. = EE a « + > ( D ° T ) - E tr {^° T ) atan ( y. ) (!3) 

! = 1 fc=l 1 = 1 ^ ! ' 

w/iere 

ai ' fe -i 1 / \ . (14) 

tr — s otherwise. 



Proof. Let us consider the expression of the T^w-norm of S, as given in Definition [2] 



: = tr (^ / l#(^)| 2 <H (15a) 



2tt 



By noting i/(s) = ^i=i ~tt¥~j it comes that 



Each term of this integral is then considered separately: 
(i) Considering the first term, it follows that 



(ii) Greater attention should be turned to the following (second) integral 

I ff(iv)n T dv> = 

2tt 



— / H{ 2 v)D T dv = — V / -^^D T dv. 



(15b) 



= tr (J- J H{jv)R{- 3 v) T dv 
As H (s) has simple poles, it can be written as 

S Ai 



l^llk. = ^ tr (/ DDT + H{]v)D T + DH{-jv) T + H{jv)H(~ ]V ) T di^j . (17) 



±- r I?# T ^ = —DD^. (18) 

27T /_,,, 7T 



(19) 



Indeed, if Re (Aj) = 0, then is integrable over [— w,w] if and only if a; < |Aj|. Hence, 

in the sequel, the following assumption is made: 

u < min (lAjl). (20) 

Re(Aj)=0 

It implies that is integrable over [— lj, lj] for all i = 1, . . . , n. Under this assumption 

and invoking formulation (39) of Definition [4] given in [Aj which states that atan(z) = 



— [log (1 + jz) — log (1 — jz)\ for z 7^ ±j, it comes that 
2j 



/ -^-D T dv = j^-[Iog(-jw-A < )-logC3w-A i )]D T (21a) 



-2^atan ( ) L> T (21b) 
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Thus 

1 f ~ . . x _ , rp 1 ^ — -S .^rp ( OJ 



/ 1 _ _ 
H{jv)D T dv = — V &£> T atan , 



(iii) In a similar way 

2 



f DH{-jv) T dv = -iVD^atan (~ 



(22) 



(23) 



(iv) Regarding the last term of (17 1, one have 



and it follows that 



1 T 1 T 

-!- / H(jv)H(-jv) T dv = YY± / — 



dv. (25) 



From here, two cases must be considered here: 
(a) If A 4 + A fc ^ 0, then 



2tt (j^ - Ai) (-ji/ - A fe ) 2nJ_ L0 \jv-\ i -jv-X k 



with p^fc = — ^ ■ Because of assumption ([20]), each term can then be integrated as 
previously and one gets 



7rAi + A fe V^i/ TrAi + Afe 
Since p^fe = p^^ the sums can be reordered as follows 



a lM = - "'7 7"' atan - . (28) 

7T Aj + Afc \\i 



(b) If A 4 + A fe = 0, then 



a ik = ± T i rr*' = ^ f (29) 



After integration, one obtains 



aa = ^f-i^V+i^x) (30a) 



1 w<fo0j[ 



(30b) 



TT w 2 + A l A l 

Taking the trace of each term and adding them up leads to the result of Theorem [T] 

□ 
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Corollary 1 (Stable and stricly proper case). Given assumptions of Theorem^ for a continuous 
stable and strictly proper MIMO LTI dynamical system of degree n with simple poles, the frequency- 
limited H2-norm is given as 



\i— 1 k— 1 



n 




tr \J2<l>iH T (-\) [ - >tan ( f ) ) ) . Cilbj 



Proof. The proof if straightforwardly obtained from ( 13 1 with Z? = and by noting that as the 
system is stable, then Aj + A& ^ 0, i = 1, . . . , n, k = 1, . . . , n. □ 



Remark 2 (About the atari term). Equation (31b I is very similar to the %2-norm expression (|5| 

presented in ' L 3, chap. 5]. Indeed, the above formula involves only the new term: — ^atan (^x~j> 

which plays the role of weight for each modal contribution of the product 4>iH T (— Aj) , as a function 
ofuj. 

Remark 3 (Imaginary eigenvalues AjJ. m ). In Theorem^ assumption is made that if H(s) has any 
imaginary poles A™, k = 1, . . . , ni m (0 < ni m < n), then uj must satisfy 

u) < min (|Ar|) . (32) 

With reference to equation (21), this assumption makes both formulations (39) and (40) in \A\ 
of the principal value of the complex arctangent, equivalent. Indeed, according to (41), the two 
formulations only differ when Re(z) = and Im(z) < — 1 which is never reached by z = 4^ due 
to hypothesis stated in relation ( 32 ) . 

Remark 4 (Meaning of the T^.w-norm for unstable systems). It is important to note that the 
TA.2.U) -norm and the H.2-norm are linked only for stable systems. Indeed for unstable systems, the 
^2,01 -norm is simply the integral (|6| and has no physical meaning. 

Proposition 1. Given a continuous and strictly proper MIMO LTI dynamical system described 
as in 0, the behaviour of the T-L2,uj-norm o/E asw->oo can be summarized as follows: 

(i) If Ti is stable, then lim Ci; _>. 00 ||E||% 2 u = ||E||^ 2 . 

(ii) If E is unstable, then lim a ,_i. 00 ||E||-H 2aj is finite (see Remark^. 

(Hi) If 12 has one (or more) purely imaginary pole(s), then lim Ll ,_>. 00 ||E||% 2 u = oo. 
Proof. The proof of each elements is given as follows, 



(i) As £ is stable and strictly proper, its %2,a;-norm is given following (31b). The limit of 
atari yx^j as w -> oo are given by considering both the definition of acot(z) and its limits 
as z -> 0, z ^ in [A] As Re(Ai) < for i = 1, . . . ,n, then 

lim atari I — J = — — for i = 1, . . . , n (33) 
u-^oo V Aj J 2 

Thus hm™ ||£|& = tr ^H T (~X l )) = ||E||^. 



(ii) Similarly as to [i), if A^ is an unstable pole of H{s) then lim^^oo atari V4A = J. By 
separating the unstable poles A^~, i — l,...,n + < n of H(s) from the stable ones A™, 
i = 1, . . . , n_ < n such that n + + n_ = n, it comes that 



Jim ||E||^ =tr ( J2<PiH T (-Xi) ) tr ( 5>tf T (-A+) ) . (34) 



,i=l / \i=l 



G 



Note that if H(s) = C(sl„ — A) 1 B has only unstable poles, then lim w _ i . 00 ||S||^ 3uj = 
ll^ _ ||%2 wnere S~ is the stable reflected system with associted tranfer function H~(s) — 
C(sl n + A)- 1 B (instead of C{sl n - A)~ X B). 



(hi) If H(s) has a purely imaginary pole Xi 
integral (191 is infinite. 



then assumption (32) is violated as uj — > oo, thus 



□ 



4 Numerical considerations &; illustration 

4.1 Numerical considerations 

The evaluation of the T^w-norm can be computed either via the Gramian-based formulation ([7| 
or via the spectral expression of Theorem [7J Given a continuous LTI dynamical system H (s) = 
C(sl — A)~ 1 B + D, the residues <fii of the system can be computed through the right and left 
eigenvectors Xi, yi £ C" of A, for i = 1, . . . , n, defined as 



Axi = XiXi , y*A = Xiy*, (35) 



using the following formula [15 



CxjgB 



, for i = 1, . . . , n. (36) 



This implies that the new spectral expression of the Hi^-norm requires the computation of 
the eigenvalues and the corresponding eigenvectors of A only. This makes the norm computation 
much easier than using the Gramian formulation. Indeed, while the latter requires evaluating 
a matrix logarithm and solving a Lyapunov equation (see Section [2]) , the spectral expression of 
Theorem [T] only requires solving an eigenvalue problem (i.e. matrix/vector operations that can 
be solved for sparse models) for which many research results provide fast and accurate procedures 

mug. 

Another advantage of the spectral expression of the Hi u -norm, once the eigenvalues and the 
residues have been computed, the T^^-norm can easily be evaluated for several values of u> with 
little extra numerical cost. This is not the case with the Gramian-based formulation where the 



last terms W c (uj) and W (uj) of the Lyapunov equations (9a I and (9b I, which have to be solved 
for each value of u, involve the computation of a matrix logarithm depending on u>. 

In order to illustrate the computational time from a qualitative point of view, the following 
test has been carried on: for a fixed u value (u> = 100), the %2,w-norm of several randomly 
generated models of dimension n = 2, . . . , 200, has been computed using each formulation^] and 
the corresponding computation time has been measured. For each model, the norm has been 
evaluated 1000 times and the mean computational time has been reported on Figure [T] as a 
function of the model order n. 

The standard deviations of the measures are negligible and have not been plotted. On this 
example, the lower complexity of the spectral formulation clearly appears in favour of the proposed 
spectral expression of the H2 iUJ -norm, especially for high-order models. 



4.2 Illustration of the behaviour of the % 2)W -norm 

Here, the Los- Angeles Hospital model provided in the COMP^ e z6 library [TU] is considered. This 
model is a SISO, stable and strictly proper model of order n = 48. Figure [2] shows the frequency 
response (top) and its frequency- limited T^-norm as a function of uj (bottom, solid blue line). Note 

2 Concerning spectral expression of Theorem [7] the eigenvalues and the eigenvectors are computed at each 
occurrence. 
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80 100 120 
Order of the model 



Figure 1: Mean computation time of the T^.u-norm for a fixed u) with respect to the order 
n = 2, . . . , 200 of the randomly generated model. Solid blue line (Spectral expression of the 
^2,u-norm, Theorem [TJ), Solid green dotted line (Gramian-based T^.w-norm). 

that here, only the upper bound of the considered frequency interval [0, u] changes. Numerical 
integration, Gramian-based formulation and the spectral formulation of Theorem [7] lead to the 
same results. This example illustrates the two following points : 

(i) As expected, for stable systems, and as uj increases, the T^^-norm increases and asymptot- 
ically tends to the %2-norm (bottom, dotted red line). 

(ii) The "H2, w -norm increases by steps. When uj crosses the abscissa of a peak in the frequency 
response (which corresponds to the absolute value of the poles of the system), the H^u-norm 
steps-up. The larger the magnitude of the peak is, the bigger the step is. This can easily 
be understood with the new spectral ^2,^-norm formula. Indeed, ( 31b[ ) is a sum of the 



product between the residues <pi and the frequency response associated to —At, weighted by 
atan yx^jt i — 1, ...,ra (see also Remark 2). It is worth being noticed that atan have a 
logarithmic behaviour between each steps. 



5 Conclusions 

In this paper, a new simple yet efficient spectral-based formulation of the frequency limited H2- 
norm, denoted %2,w-norm, has been presented. For stable and strictly proper systems, this spectral 
formulation is very similar to the %2-norm presented in [3] . More especially the newly introduced 
arctangent term plays the role of weights on the contribution of each modal content and gen- 
eralizes this formula. In addition to be an alternative way to formulate problems involving the 
%,w-norm (such as controller design, robust analysis or in model approximation...), this spectral 
expression also offers some numerical advantages over the Gramian-based formulation ([7]). Indeed 
while the latter requires evaluating the matrix logarithm and solving a Lyapunov equation, the 
former one only requires the eigenvalues and the eigenvectors computation, which is numerically 
less expensive. Moreover, this new formulation makes possible the computation of the ^ 2 .u- 
norm to large-scale systems insofar as sparsity is exploitable and there exists efficient iterative 
algorithms to get eigenvalues and eigenvectors (see e.g. [TH [T3]). To the best of the authors' 
knowledge, this spectral formulation may help reconsidering some control-oriented problems from 
a new perspective and provide a consistent tool for many research and applicative topics. 
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Bode Diagram 




of the frequency interval 



Figure 2: Frequency response of the Los- Angeles Hospital model (top) and its %2,w-norm (bottom) 
as a function of w. 



A Complex function definition and conventions 

This appendix is mainly based on [5] and is dedicated to the presentation of the definition and the 
conventions which has been chosen for some complex functions involved in the proof of Theorem 
[T] Given a complex number z, the principal value of its argument, denoted as arg(z), is defined 
such that 

— 7r < arg(z) < 7r. (37) 

The function arg(z) is continuous everywhere in the complex plane excepted on a branch cut 
along the negative real axis. Each function listed below inherit its branch cuts from the principal 
value of the argument. 

A.l Complex logarithm 

Definition 3 (Principal logarithm). The principal value of the logarithm of z, denoted log(z) ; is 
defined for z =/= as 

log(z) =ln(|z|)+jarg(z), (38) 
where ln(x) is the natural logarithm of x € W + . 



A. 2 Complex arctangent 

Definition 4 (Arctangent). The principal value of the arctangent functions can be defined in two 
different ways denoted here atani(z) and atan2(z). For z ^ ±j. 

atamO*) = — [log (1 + jz) - log (1 - jz)\ (39) 
2j 



and 

atariof z) = , 

2j Vl-.^ 



atan 2 (z) = ^log ( \^ ) . (40) 
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It can be shown that atani(z) = atan2(z) everywhere in the complex plane excepted on the 
branch [—j,—joo) [5], indeed 

I ^ (l+jz\ J 7r+^[log(l+.?z)-log(l-jz)] if Re(z) = and Im(z) < -1 



2 3 V 1 3* J | _L [log (1 + jz) - log (1 - jz)} otherwise. 

(41) 

A. 3 Complex arccotangent 

Definition 5 (Arccotangent). The definition of the principal value of the arccotangent of z, noted 
acot(z), is based on the definition of the arctangent function, indeed for z ^= ±j and z =/= 0, 

acot(z) = atan (~J- (42) 



Definition [5] leads to two formulations for the arccotangent function, acoti(z) = atani(^) and 
acot2(z) = atan2(i), depending on which formulation of the arctangent is considered. These two 
formulations are equivalent excepted on ]0,j[. Indeed for z ^ ±j and z =/= 0, 



2; 



log 



z + J 
z-j 



7T H 

1 

27 



1 

27 

log! 1 



log 1 



logf 1 

J 



- log 1 - 



if Re(z) = and < Ira(z) < -1 
otherwise. 



Note that acot2(z) is defined for z — whereas the limit of acoti(z) when z 
Hence 

f for Re(» > 



hm acoti(z) 

2— >o 



f for Re(z) > and Im(z) < 
■| for Re(z) < 
■f for Re(z) < and Im(z) > 



(43) 

is not unique. 



(44) 



which can be simplified due to equivalence ( 43 1 in 

lim acoti(^r) = 



z->0 



for Re(z) > 
for Re(z) < 0. 



(45) 



In the proof of the main result (Theorem [lj, formulation (39 1 is used to define the complex 
arctangent, but, due to theorem hypothesis, formulation (40) would have led to the same result. 
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